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Abstract 

In bricklayers' model, which is a generalization of the misan- 
thrope processes, we show that a nontrivial class of product distri- 
butions is closed under the time-evolution of the process. This class 
also includes measures fitting to shock data of the limiting PDE. In 
particular, we show that shocks of this type with discontinuity of 
size one perform ordinary nearest neighbor random walks only in- 
teracting, in an attractive way, via their jump rates. Our results are 
related to those of Belitsky and Schiitz [4] on the simple exclusion 
process, although we do not use quantum formalism as they do. The 
structures we find are described from a fixed position. Similar ones 
were found in Balazs [2], as seen from the random position of the 
second class particle. 
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Introduction 

It has been known for a while that many models from the field of stochastic 
interacting particle systems have deterministic partial differential equations 
as scaling limit, see e.g. the book of Kipnis and Landim [11], or Seppalaincn 
[13], Toth and Valko [15]. In the so-called Eulerian scaling this equation is 
usually a Burger-type one, developing shocks as entropy solutions from any 
decreasing initial data. The microscopic structure of these shocks is of great 
interest. Many results construct such structures from the viewpoint of the 
so-called second class particle, which is an object coming from coupling of 
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these models. Some examples are De Masi, Kipnis, Presutti and Saada 
[6], Derrida, Lebowitz and Speer [7], Ferrari [8], and Ferrari, Fontes and 
Kohayakawa [9]. As this particle has a complex random motion in the 
system, it is not immediately clear how these structures look like from a 
fixed, non-moving position. On the other hand, this way of looking at the 
problem makes it difficult to deal with the case of multiple shocks. Ferrari, 
Fontes and Vares develop these methods and handle this case in [10]. 

Another work in this area is by Balazs [2], who introduces a simple 
product measure which shows the properties of a shock, and is stationary 
as seen by the second class particle. This result is valid for the exponential 
bricklayers ' process (see below for a definition) . Our aim here is to identify 
a structure corresponding to this distribution, but this time as seen from 
a fixed, non moving position. This way of looking at the problem makes 
it easy to handle multiple shocks as well, thus we can also investigate how 
these shocks interact with each other at the microscopic level. 

We consider the class of bricklayers ' processes, introduced in [2] and 
[3] based on ideas of B. Toth. These models are slight generalizations of 
the misanthrope (see Cocozza-Thivent [5]) and also of the zero range pro- 
cesses. Each model from this class has a special one-parametered family of 
one dimensional discrete measures. If we fix that parameter and build the 
product of these measures for the sites, then we obtain a time-stationary 
distribution of the model. Hence we have a one-paramctcrcd family of equi- 
librium distributions; this parameter sets the slope of the wall, a quantity 
corresponding to particle density in particle systems. 

We consider product measures with the same marginals as for the sta- 
tionary distribution, except for that we allow the parameter to change from 
site to site. The main result states that distributions in this class evolve 
to linear combinations of distributions from the same class, but only in the 
special case of exponential bricklayers' process. The form of the linear co- 
efficients allows us to interpret some of the situations as ordinary random 
walks of shocks having discontinuity of size one. We also obtain the nature 
of interaction between these shocks. It follows that a group formed by a 
number of such one-sized jump shocks is of stochastically bounded size in 
time, i.e. shocks of larger (integer-valued) discontinuities represented by 
such a group are sharp under any kind of hydrodynamic scaling. 

Bclitsky and Schiitz [4] derive results similar to ours for the simple exclu- 
sion process with the use of a quantum algebra symmetry Unfortunately, 
we were not able to develop the quantum formalism for our locally infi- 
nite state space systems. However, there is a remarkable analogue between 
their work and the present settings: in [4] special relations are required 
between the shock densities and the particle jump rates. These relations 
are identical to those of Derrida, Lebowitz and Speer [7] allowing there an 
exact (not only asymptotic) product-description of the stationary distribu- 
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tion as seen from the second class particle. In the present work we shall 
be specially requiring integer values for the size of discontinuities of our 
shocks, in a very similar way as needed for the exact result of Balazs [2] 
from the viewpoint of the second class particle. 

The model and its stationary distributions are introduced in the next 
section. Then we give a very brief overview on the Eulerian hydrodynamic 
limit of these models in Section 3. Our results are contained in Section 4, 
and the proofs can be found in the last section. 

2. The bricklayers' model 
2.1 Infinitesimal generator 

In [3], based on the idea of B. Toth, Balazs defines a class of deposition 
models including the well-known simple exclusion and zero range processes. 
The so-called bricklayers' model to be discussed below is a family from the 
class described there. We briefly repeat the definition here, see [3] for a 
more detailed introduction. 
We consider the state space 

(1) n = {u = (wj)iez : w ; eZ} = Z z . 

For each pair of neighboring sites i and i + 1 of Z, we can imagine a column 
built of bricks, above the edge (i, i+1). The height of this column is denoted 
by hi. If co(t) € f2 for a fixed time t £ K. then u>i(t) = hi-\(t) — hi(t) e Z is 
the negative discrete gradient of the height of the "wall" . We are basically 
interested in this quantity, not the absolute heights h of columns. Given 
a configuration, the growth of a column is described by a Poisson process. 
A brick can be added to a column, with rates depending on the local 
configuration: 

lu — > o/ hl+1 ^ with rate r{uJi) + r(— Wj+i), 

where 

=(•••, Ui - 1, Ui+l + 1, UJ l+2 , ...) 

only differs from u at sites i and i + See fig. 1 for some possible instan- 
taneous changes. The process can be represented by bricklayers standing 
at each site i, laying a brick on the column on their right with rate r(w,) 
and laying a brick to their left with rate r(—Ui). This interpretation gives 
reason to call these model bricklayers' model. 

The rate function r is a positive real valued function on Z. In order to 
have product measures as stationary ones (see Corollary 5.2), we assume 
that r satisfies 

(2) r{z) ■ r{-z + 1) = 1 for all z e Z. 
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r(u)i) + r(-w i+ i) 



Figure 1: A possible move 
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At time t, the interface mentioned before is described by uAt). Let 
ip : Q, — > R be a bounded cylinder function i.e. ip depends on a finite 
number of variables Wj. The growth of this interface is a Markov process, 
with the formal infinitesimal generator L: 

(3) {L v )(u) = £{[r( Wi ) + r(-w i+ i)] ■ - <^)] }• 

Note that for each index i, uji can also be negative hence direct particle 
interpretation fails. Apart from this, L is formally the sum of two totally 
asymmetric zero range-generators; one for u> as number of particles which 
jump to the right, and another for — ui as number of particles with jumps 
to the left. Therefore, part of our arguments can easily be applied to zero 
range. We shall notify the reader where the analogy ends. Unfortunately, 
the most interesting conclusions do not apply to the zero range process. 

When constructing the process rigorously, problems may arise due to 
the unbounded growth rates. Existing methods for constructing such pro- 
cesses only apply when the rate function r does not grow faster than linear, 
see Andjel [1], Booth and Quant [12]. However, models with exponential 
rates r are of essential interest for this note. In a forthcoming paper we 
prove that, in the attractive case, when r is a monotone increasing func- 
tion and does not grow faster than exponentially, the growth of any fixed 
column until any fixed moment is a random variable of finite second mo- 
ment if we start the process from an appropriate set of initial configura- 
tions. Unfortunately we are not able to carry out a full construction of this 
Markov-process at this moment. We assume thus that existence of dynam- 
ics on a set of tempered configurations ft (i.e. configurations obeying some 
restrictive growth conditions) can be established. Technically we assume 
that il is of full measure w.r.t. the product measures defined in Section 2.2. 
We do not deal with this question in the present paper. 

The exponential bricklayers' model 

A special case of the models is the exponential bricklayers' model (EBL), 
where for z G Z 

(4) r(z) = e-£ e? z 
with a positive real parameter /?. 

2.2 The product measure 

For n G Z + , we define 

n 

r(n)\ := JJ r(y) 
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with the convention that the empty product has value one. For general 
z e Z, we then define 

r{z)\ : = r(|«|)!. 

Using (2), the property 

r(z + 1)! = r{z)\ ■ r(z + 1) 

can be shown for any integer z, giving justification of the notation used 
here. Let 

9 : = liminf log(r(n)) 

n — >oo 

which we assume to be strictly positive. By (2), it is in fact strictly positive 
in case of monotonicity of r or, in other words, attractivity of the model. 9 
can even be infinite. With a generic real parameter 9 e (— 9, 9) , we define 



Z{6):= J2 



and the measure 

(5) ^(z):- ' ' 



Z{9) r{z)\ 

on Z. The definition leads to two simple but essential properties: 

^(z-1) 

m » (6){z) 

6 ' WW) =r( ' z) ' 

of which the latter is a consequence of (2). Based on these, we also see 

(7) Er(z) = e 9 and Er(-z) = e" fi 

w.r.t. 

For the special case of the EBL model, for 9 E (— oo, oo), we obtain the 
discrete normal distribution 

8 ^ (8) W - — = ~,_ , 

e -WZ(9) Z{[3, m) 

with the notation m := 6/(3. 

For the purposes of this paper, we introduce the parameter-vector 

0:=(0i)i€Z where 0; e (-9, 9), 

and the product distribution on f2, having marginal /x^-' for site i. We 
also assume that 9i has finitely many different values as i goes from — oo 
to oo. 
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3. Hydrodynamic limit 



As we shall see later, the product measure jjf^ is stationary for the models 
if Oi has the same value for all sites i. Allowing this equilibrium to change on 
a macroscopic level, the hydrodynamic limit of these models can be derived. 
For different types of scaling limits we refer to Toth and Valko [15]. Here 
we only consider the so-called Euler-scaling. Strictly speaking, as in case of 
construction of these processes, we only know about the case of sublincar 
r functions handled in the context of hydrodynamic limits. However, for 
this note exponential growth rates are essential. Via formal computations, 
one can guess for the form of the Eulerian limiting hydrodynamic partial 
differential equation even in this superlinear case. The variable of this 
equation is formally defined by 

(9) ti(i,i): = Eu l/e (t/£). 

Then via formal computations we obtain the differential equation 

, „. du dJ(u) 

as e — > 0, where J(u) is defined as follows. The function u{6) = (uj) of 
9 is strictly increasing since the derivative 



du(fl) 
dO 



E 



(*V)- (eW(w)) : 



is positive (—6 < 9 < 6). Let 9(u) be the inverse function. The quantity 
E( 9 ) (r(uj) + r{—oj)) depends on 9, and J is defined by 

(11) J{u) := EW M » (r(w) + r(-w)) = 2 cosh(6»(w)). 

In [3] , Balazs proves that this function is strictly convex if r is monotone 
increasing and is convex. (An example of such rates is the exponential 
bricklayers' process defined above.) In this case, (10) develops shock (weak) 
solutions. A (single) shock solution starts with initial data 



ftileft , X < 
[bright , X > 



u(0, x) 

with ui e ft > u r ight- The stable weak solution is of the form 

u(t, x) = 



fltleft , X < St 
\ Wright , X > St 
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where the speed s of the traveling shock is determined by the Rankine- 
Hugoniot formula 

( 12 ) g = Jjuicit) - JKight) 

Mloft - Wright 

see e.g. Smoller [14]. We define the size of the shock as ui e ft — Wright- 

In the case of multiple (but finitely many) shocks, i.e. piecewise con- 
stant, decreasing initial data with more than one discontinuities, a similar 
formula holds for each shock if u has values ui e ft and w r i g ht at its two sides. 
By convexity of J, these shocks meet in finite time, forming less and less 
shocks of larger and larger sizes. 

Using definitions (11), (9) and the fact (7) in the Rankinc-Hugoniot 
formula (12), the speed of a shock can now be written as 

p0(ltleft) I p-0(«left) _ p0("right) _ p-S(«right) 

(13) 8 = - — . 

Wleft — Wright 

This is what we see on the macroscopic scale. Our aim now is to find 
some similar structures in the microscopic level. 



4. Results 

For simplicity, we introduce the notations 

^ '■=(■■■, 0» + 0, ■ • ■ ) 

8^' ^ :=(..., di-i, $i — /3, i+ i, ... ) 

given a parameter (3. 

Theorem 4.1. For the exponential and no other bricklayers' model, (nei- 
ther for any totally asymmetric nearest-neighbor zero range process), the 
product measure defined at the very end of Section 2, evolves to a lin- 
ear combination of product measures with similar structure. The evolution 
is formally described by the equation 

( 15 ) ^H m = E (e e --e^)-(^ (l+1 ' +,) -M®) + 

+ E (e-^-e-^).(^ (i - 1 '- , )-^ ) ). 

The proof can be found in Section 5. We cannot give a nice interpre- 
tation of this evolution in the general setting. However, in some special 
cases, this theorem gives rise to a rather interesting observation. For the 
rest of this section, we are talking about the case of exponential bricklayers' 
process. 
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4.1 One single (and special) shock 

Fix a parameter loft for the left-hand side, and # r i g ht : = #i e ft ~ P- Fix a 
position Q eZ and 



9,: = 



for i < Q - 1, 
?ri g ht for i > Q 



for components of 9. Then is a shock-measure corresponding to a shock 
of size one: for the exponential model, a change of size (3 in the parameter 9 
implies a jump of size one in the expectation u = E^^(w) (see also Lemma 
5.3). We say that this shock-measure is located at position Q — \. Balazs 
[2] shows that this measure is in fact stationary as seen from the random 
position of the so-called second class particle. In the present setting, we 
examine this structure from a fixed position. First, we apply (15) to this 
Cfi). 



(16) A.f/0) = ( e *e ft _ e W») . (> <Q ' +) ) - + 

+ (e-W._ e -^).(/ £ (fi w - 1, - ) )- H (e). 

Then we realize that 

q(Q,+) _ f fteft for * < Q, 

1 \ bright for i > Q + 1, 

and 

/,(Q-1, -) _ / fteft for i < Q - 2, 



bright for i > Q - 1. 

The former is the same shock measure located at <5 + |, the latter is located 
at Q— 1|, see figures 2 and 3. Thus (16) says the following: a shock measure 
of size one, located at Q — | becomes a similar measure at Q + | with rate 
c #icft _ e fri g ht ; or becomes one at Q — 1| with rate e _e "s ht — e _6,loft . Hence the 
shock performs an ordinary nearest-neighbor random walk with constant 
rates e e ' cft — e e " sht to the left and e~ 6risht — e~ eic!t to the right, respectively. 
These rates are positive according to 0i e ft > r ight) and they in fact agree 
with the expectation of the left and right jump rates of the second class 
particle w.r.t. fi^-\ The expected speed of this random walk therefore 
agrees with the speed of the macroscopic shock (of size one) suggested by 
the Rankine-Hugoniot formula (13). This is the shock-measure of [2] as 
seen from outside. 

We emphasize that the ordinary random walk shown here and the sec- 
ond class particle of [2] are different objects. While the jump rates of the 
latter depend on the actual state of the model, in the present context we 
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> 



e r,gh 



9i 

"left " 
;hfT 



9rl 



H h 



Q-i Q 



H 1 h- 



H 1 1 1 ► 1 

the random walker ^+2 



H 1 1 1 ± 

the random walker ^+2 



Figure 2: Step of the random walker from Q — ^ to Q + \ 



01eft " 
bright" 



> 



•ight _ p — #left 



H h 



0< 

#lcft " 
bright" 



Q-i Q 



> 



H 1 h- 



H 1 1 1 1 

the random walker z +2 



•«■-- 

H 1 1 1 «- 1 

the random walker z +2 



Figure 3: Step of the random walker from Q — \ to Q — 1^ 
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only deal with distributions of the process, hence we can not even state any 
kind of dependence of our random walks on particular states of the model. 
One could possibly think about comparing the law of the second class par- 
ticle, integrated out w.r.t. the initial distribution fi^S) of the model, to the 
distribution of our random walkers. We do not know about any argument 
which could bring these laws in connection to each other. In particular, 
by averaging out the continuously changing random environment this way, 
nothing assures a priori the Markov-property of the second class particle, 
while our ordinary random walkers clearly have this behavior. 



4.2 Multiple (special) shocks 

As for coordinates of 0, assume that 0j as a function of i 

• is decreasing, 

• changes finitely many times 

• in multiples of (3. 

We also assume that the values of 0j range from r ; g ht to 0i e ft = r i g ht + n/3. 
As each jump of size (3 in the parameter results in a jump of size one in 
the expectation (see Lemma 5.3), and so in the hydrodynamic quantity u, 
now corresponds to a multiple shock measure, where each shock has 
integer size giving a sum of n. In a similar way as for the single-shock case, 
one reads the following representation of (15). If 0r = 0l — m/3, and 



(17) 



9 L fori = Q-l, 
9 R for i = Q, 



then we have a shock of size to located at Q — |. With rate e SL — e 6n , 
0q increases by (3, i.e. the size of the shock at Q — \ decreases by one and 
the size of the shock at Q + \ increases by one (or a new shock of size one 
emerges if there wasn't any shock here before). With rate e~ 0K — e~ 6lj , 
0q_i decreases by (3, i.e. the size of the shock at Q — | decreases by one 
and the size of the shock at Q — \\ increases by one (or a new shock of 
size one emerges here). See figures 4 and 5 for these steps. This, in fact, 
gives us the following picture: the n one-sized shocks can be represented 
by n random walkers, each performing a nearest-neighbor random walk on 
Z + | . These walkers meet eventually, hence any number to < n of them 
can be present at a given position. Assume we have to > walkers at this 
given position, and k of these walkers strictly left to this position. In this 
case, since each walker (i.e. each shock) changes the value of by (3, we 
have 0l = 0ieft — k(3 and 0r = 0i e ft — (fc + m)(3 with the notation of (17). 
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H h 



<3-i Q 



H 1 h- 

Q-l Q 



H 1 1 1 ►- 1 

random walkers z +5 



H 1 1 1 ►- 

random walkers z " 



Figure 4: A shock of size one steps to the right from position Q 
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> 



>2/3 



e -0R _ c -9l 



H h 



Q-i Q 



H 1 1 1 >- 1 

random walkers z +5 



>2/3 



> 



— I 1 1 — 

Q-i Q 

•<--- 

H 1 1 1 

random walkers z + 



Figure 5: A shock of size one steps to the left from position Q — h 
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Therefore, the rate with which any walker steps from here is 

e 0L _ e 0R = e e lott -kl3 . I X _ e -m0\ t the ri ht and 

18 v ; 

e -e R _ e -e L =e -e left +fe/3 . ( e m/3 _ ^ t0 the left) 

respectively. If a walker is alone at her site and she has I > walkers to 
her left, then we have m = 1 and write I instead of k. In this case, our 
formulas reduce to e eieH ~ l @ (l — c~^) as for rate to jump to the right with, 
and e~ e ieH+W ( c /3 _ ^ as f or ra ^ e ^ j um p to the left. In fact, the jump 
rates (18) are just the sums of the latter jump rates for the single walkers 
as I goes from k to m + k — 1. In other words, the rate with which any 
walker jumps from a position is just the sum of those rates if the walkers 
at that position were alone. 

The expressions for these rates show that walkers on the left (smaller k 
values) have larger rates to jump to the right and smaller rates to jump to 
the left than others more on the right (larger k values). Only their order 
matters, not their absolute position. The form of these jump rates is the 
only way the shocks interact. 

We also see from here that the shock size at any fixed position can only 
change by one at a time. Hence shocks of size larger than one can never 
move, only shocks of size one are separated and merged again eventually 
somewhere else. On the microscopic level, larger shocks are only formed by 
accidental meeting of shocks of size one which perform the random walks. 

The rate of these walkers (or shocks) agree with the expected jump rates 
of the second class particle w.r.t. ^-K However, as mentioned above for 
the case of a single shock, the second class particle and our random walkers 
are conceptually different objects. 

As larger shocks can only move by separation and displacement of one- 
sized shocks, one can ask if these several little shocks in fact show the 
properties of one large shock on the macroscopic scale. An answer to this 
question is given below: the following lemma shows that these type of 
shocks are sharp under any kind of hydrodynamic limit. We define the 
width of such a larger shock (formed by a number of one-sized shocks) as 
the difference between the position of the leftmost and the rightmost shock 
(or random walker) . 

Lemma 4.2. The width of our types of shocks has bounded moments in 
time. 

The proof of this lemma is contained in Section 5. 

As the one-sized shocks (i.e. walkers) stay in stochastically bounded 
distance from each other, a natural quantity indicating the position of this 
group of walkers is their central mass. To be more precise, let X^(t) 
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denote the position of the m-th walker at time t, m = 1 . . . n (with any 
labeling), and define 

(19) X(t) : = X W (t)+XW(t) + --- + XW(t) _ 

Lemma 4.3. The process X(t) is an ordinary continuous-time random 
walk on (Z + 1/2) /n with 

right jump rate e e ' cft — e SriEht , 
left jump rate e ~ e "« ht - e~ 01clt . 

The proof of this lemma can also be found in Section 5. As a corollary, 
X(t) satisfies strong law of large numbers 

X (t) e^ lott + e~^ loft — e^ right — e — " r 's ht 
lim = a.s. 

t^oo t n 

As the n walkers together represent a shock of jump size n, this formula 
is in accordance with the Rankine-Hugoniot formula (13), derived for the 
speed of shock solutions of the PDE. 

The model and the situations detailed above are the ones for which we 
could establish existence of these kind of solutions in terms of ordinary 
random walks. Situations when the random walk interpretation fails are 
listed below. 



4.3 Other initial conditions 

Unfortunately, we don't have such a nice interpretation when any of the 
conditions in Section 4.2 fails. Theorem 4.1 still holds in this case as it 
does not depend on these conditions, hence the product structure of the 
measure is still closed under the time-evolution of the system. One just 
cannot give the random walk representation of the result. If 6i as function 
of i is not decreasing at some i value, then the rates of our random walkers 
would become negative; we don't know how to interpret this case. It is not 
surprising that such random walk description fails: upward jumps of the 
function 6i correspond to upward shock solutions. By strict convexity of 
the hydrodynamic flux function J(u), which is in fact established in [3] for 
the exponential rate functions r we have in the exponential model, upward 
shocks are non stable, they evolve into rarefaction waves (see e.g. [14]). 
Hence it is not reasonable to expect some kind of random walks of upward 
jump solutions. 

On the other hand, if 6i is decreasing but not in multiples of /?, that 
means that we have shocks of non- integer size (see Lemma 5.3). By sep- 
aration of one-sized shocks, as described by the Theorem, eventually a 
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shock of size smaller than one appears. When another shock of size one 
jumps from its position, it leaves an upward shock behind itself. At this 
moment, our walkers' rates become negative, and the interpretation fails 
again. Unfortunately we don't have an intuitive argument explaining this 
phenomenon. 

4.4 The not completely asymmetric case 

Bricklayers' process in the presented form corresponds to totally asymmet- 
ric particle systems. One can modify the generator by multiplying the rates 
so far introduced with a factor p and include brick-removals as well: 

lo — with rate q (r(— Wj) + r{uJi + \)) , 

where 

W (l+M) =(••., Wi + 1, UJ i+ l - 1, L0 l+2 , ...) 

and q+p = 1. Theorem 4.1 can easily be extended to this case: the class of 
product measures with marginals /j,^ of spatially changing 9 parameters 
is closed under the time-evolution of the process. The question is if the 
random walk interpretation can also be applied to this more general case. 
While well-posed asymmetry of simple exclusion is essential in [4] , it proves 
to be useless for our random walk interpretation in the bricklayers' process. 
Independently of the initial distribution (cither increasing, or decreasing), 
general asymmetry of the model would immediately imply negative jump 
rates of our random walkers when trying to introduce them. Hence we 
cannot generalize our random walk interpretation, or weaken any of the 
conditions of Section 4.2, needed for this interpretation, by introducing 
general asymmetry in the system. 

4.5 Other models 

As seen above, other bricklayers' models do not even show closedness of the 
product measures under the time-evolution as exponential bricklayers' pro- 
cess does. For general bricklayers' models or shock data, we do not expect 
qualitatively different behavior, rather results similar to the exact ones ob- 
tained here to hold asymptotically in space. For example, based on results 
obtained in this direction for simple exclusion (see e.g. Derrida, Lebowitz 
and Speer [7] or Ferrari, Fontes and Kohayakawa [9]), we conjecture that 
the stationary distribution as seen from the single one second class particle 
converges asymptotically to the product measure fj,^> with the appropriate 
#i e ft and bright values. We do not have any results in this direction. There 
is only an easy argument in Section 4.2 of Balazs [2] showing that once this 
convergence and also ergodicity of the stationary measure is established, 
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the second class particle satisfies the Law of Large Numbers with the speed 
value given by the Rankinc-Hugoniot condition. 



5. Proofs 



For simplicity, we introduce the notations 



(20) 



;(*,+) 



(. . . , LJi-1, LUt + 1, U)i+\, . . . ) 
(. . . , LUi - 1, LJi+1, . . . ). 



Lemma 5.1. Let ip be a bounded cylinder function on £1. Then according 
to the bricklayers' dynamics with any rate function r (satisfying (2)), 



(21) 



dt 



. E GD 



vGsO - E - e9,+1 ) ' ( E(iD ^ (<+1, +) ) - E( - } <fte)) + 



+ Y, (e- fl '-e- fl '- 1 )-(E^v'(w (i " 1 '" ) )-E^^(w)) 



where the expectation E® is taken w.r.t. the product measure /i® defined 
at the end of Section 2. 

Proof. Using the generator (3), we start computing this time derivative. 
By (6) we get to 



E® ip(w) = E^- J (Lip){w) 



dt 



= E® J ^ ([ r ( w< ) + r(-w i+1 )] • [^ (M+1) ) - ^] ) } = 
Uez J 

x <^ (M+1) ) -<p(w) 



Under this expectation, the fractions above increase u>i and decrease Wj+i in 
the argument of ip by one, respectively. (This can be seen e.g. by changing 
those variables under E®.) Remembering what uA l ' l+1 \ u/- z ' + \ uA l '~^ 



1G 



means, this leads us to 



^ V (o0 =E^){X: [e fl V(w (i+1 ' +) ) - e e V(w (< ' +) ) + 

+ e" ei +V(w (i ' ~ } ) - e -°*+V(w (i+1 ' _) ) 

Now let I and r denote the indices for which the cylinder function <p depends 
on only wj, . . . , w r . Then ip(u/- % ' + ^) and </j(o/*' ~)) are identical to for 
i's not between Z and r inclusive. Hence 

-^E^^) = E®{ ^ fe fl V(w (<+1 ' +) )-e fl V(w (<,+) )+ 

dt U=;-i 

+ e^'+V^' - } ) - e - 9i+ V(w (<+1 ' _) )] | = 
U=;-i 

+ [e-**- e -*-]^-i>-)) + 

+ e e ^(w r+1 '+)-e 9 '-V(w i - 1 ' + ) + e- e V(w i - 1 '-)-e-^+V(w r+1 ' + )| = 

= E<®{ ^ [e e * - e^ 1 ] 
U=i-1 

+ ^ [e-^-e-^]^- 1 '-))+ 



i=J+l 



+ [e 9r -e"'" 1 +e" 



] 



+ 



r+1 



i=2+l 



This formula only differs from the right-hand side of (21) in indices i where 
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the summands are zero (either because 0's agree or due to the cylinder 
property of ip). □ 

Corollary 5.2. The product measure with 6i = constant is time- 
stationary for the bricklayers' process. 

Remark. Note that up to this point a computation for the totally asym- 
metric nearest-neighbor zero range process would look very similar. For 
these arguments to apply to zero range, one would have begun with uji G Z + 
in (1), and put zero to the right-hand side of (2). Then formally we obtain 
our formulas for zero range by simply neglecting all terms containing r(— ui) 
and c~ e . However, the next lemma requires doubly infinite possibilities for 
values of u>i, hence, by its essence, it cannot be applied to zero range. 
Fix 9 and 

(22) 6' :=(..., Bi-u 0', i+ i, ...)• 

We are after the possible values of 9', if any, for which E^ip(u/ 1 ' ± ^) can be 
written as the expectation of y(u>) w.r.t. the expectation E^- ) according 
to this new parameter- vector 0'. 

Lemma 5.3. With the notation of (22), 

(23) E^V(^ (4,±) ) =E^'V(o;) 

for any cylinder function ip if and only if the model is the exponential 
bricklayers' model with parameter (3 — ±(8' — 9i) (see (4)). In other words, 

= E<* (<, %( £! i) 

holds only for the exponential bricklayers' process, and there is no other 
set of parameters which can be written in place of 9^'^ (see (14) for the 
definition of 9 {t ' ±] ). 

Proof. As E® and E^- ' are product measures, (23) is equivalent to 

^'(^1) = /'W 

for all z G Z. By the definition (5) of these measures and property (2), we 
can expand this equation: 

I e «,(zTl) i e 6'z 



(24) Z{9i) r(z T iy. Z(9>) r(z)\ 

r (± z ) = Z ^ . e (.8'-Qi> e ±e « 



18 



From this we see that r cannot be other than exponential of the form (4) 
with parameter (3 = ±(8' — 9i). We need to verify that in this case (24) in 
fact holds. To this order, we assume 9' — Oi ± (3 and r(z) = e _ ^/ 2 + /3z . We 
observe that due to the discrete normal structure, Z(f3, m) defined in (8) 
is periodic in the variable m with period one. This means 

= e T»'+4e«f Z (ft = <*'+i Z(fl') = Z(fl'). 

Plugging this into the right-hand side of (24) we have 

. Jfi'-O^z ±0t _ -g (fl'-fl 4 )* _ e -f +/3(±z) 
Z(6>) c -e e 

which agrees to the left-hand side of that equation. □ 

With this result and definitions (14), (20) in mind, (21) can be written 
in the form 

^E® = ( e * - ' ( E( - +>) <Pl& - E( " } ^) + 

+ E - • ( E( -'"' _>) <p&) - E( - } ^)) 

for the exponential bricklayers' process, and no other bricklayers' (nor to- 
tally asymmetric nearest-neighbor zero range) processes. Hence we have 
arrived to the proof of Theorem 4.1. 

Remark. For other than exponential rates, Lemma 5.1 tells us that jjS-^ 
evolves into a product of our original marginals ii e {-) and of these marginals 
shifted upwards or downwards by one: fi e (- ± 1). However, these shifted 
measures do not correspond to the original marginals with shifted param- 
eters as Lemma 5.3 does not apply for other than exponential models. 
Therefore the class of product measures of marginals [i e is not closed un- 
der the time-evolution in this case. As we do not know how the product 
of shifted versions (in the sense above) of // evolves, we do not have any 
positive result for the non-exponential case. 

Proof of lemma 4-2. For the purposes of this proof, we enumerate all n 
random walkers. At each moment, we label them from 1 to n in such a 
way that a walker right to another must have larger label than the other 
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one. We translate their jump rates according to this labeling: while the 
ra-th walker is alone at his site, she has rate e e ieft-(m,-i)/3 _ e e lo{t -mf] to 
jump to the right, and rate e" (eioft ~ m/3) - e -(.^eft-(m-i)0) to j ump to tlle 
left, respectively. If more of them meet at a site then the lowest labeled 
of them has the sum of these left-rates as rate to jump left with, and the 
largest labeled of them has the sum of the right-rates to jump to the right 
with, respectively (see the remark after (18)). All others have zero rate to 
jump due to our labeling method. These rates arc in accordance with the 
rates defined in Section 4.2. 

Define rj m as the distance between the m-th and the m + 1-th walker 
(m = 1 . . .n — 1). Then r/ m > 0, and, whenever it is positive, it decreases 
by one with rate larger than or equal to 

e e lott -(m-l)/3 _ gSieft-m/3 + e -(e lcft -(m+l)/3) _ e _(0 left - m /3) _ 

= c' 3 (1 - c-P) (e e >^ m + e - e ^ +l3m ) . 

This would be the rate if the m-th and m + 1-th walkers were alone at their 
sites. If they are not alone, then by the labeling method, to is the first one 
to jump to the right and m + 1 is the first to jump to the left, hence rj m 
decreases even faster. However, when r\ m is positive, its maximum rate to 
increase is 

e -(0 left -m/3) _ e -(0ieft-(m-l)/3) + e loft -m/3 _ e left -(m+l)/3 _ 

= (l - e' f3 ) ( e eioft ~' 3m + e~ eirft + /3m ) . 

This would be rate for i] m to grow if the walkers to and m + 1 were alone 
at their sites. If they arc not, then by the labeling method to is the last 
one to jump left and to + 1 is the last one to jump right, so this rate is even 
smaller. 

Thus we see that i] m is a process on Z+ which increases with smaller 
rates than decreases (except when it is zero) . Hence its stationary distribu- 
tion is stochastically bounded by a geometric one (disturbed at zero) and 
so all its moments are finite. As the width of the shock is just 

n-l 
m—l 

we see that it also has finite moments. □ 

Proof of lemma 4-3- Jumps of walkers from different positions happen in- 
dependently, with rates 

c e L _ 6 0r to- the right, and 
e -e R _ e -0 L t0 tne left 
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with #l and #r being the local parameter values to the left and to the 
right of the position in question. The rate with which any of our n walkers 
jumps is therefore the sum of the previous rates, and is 

, . e 6, " il — e e "z ht to the right and 

' ' e -eright _ e -flieft to the left, 

where #i e ft (Aright) is the value of the parameter left to the leftmost (right 
to the rightmost, respectively) walker. 

X(t) defined in (19) is a rescaled sum of the walkers' position, therefore 
it increases (decreases) by 1/n if and only if any of our walkers jumps to 
the right (left, respectively). This happens according to a Poisson-process, 
with rates shown in (25). □ 
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